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Abstract 



We show how to define and calculate the ground state energy of a sys- 
tem of quantum particles with 5 attractive interactions when the number 
of particles n is non-integer. The question is relevant to obtain the prob- 
ability distribution of the free energy of a directed polymer in a random 
medium. When one expands the ground state energy in powers of the 
interaction, all the coefficients of the perturbation series are polynomials 
in n, allowing to define the perturbation theory for non-integer n. We de- 
velop a procedure to calculate all the cumulants of the free energy of the 
directed polymer and we give explicit, although complicated, expressions 
of the first three cumulants. 

PACS numbers: 05.40. -fj, 02.50.-r, 82.20.-w. 

It is our great pleasure to dedicate this work to our friend Joel L. Lebowitz on 
the occasion of his 70th birthday. 

1 Introduction 

We consider a system of n identical quantum particles on a ring of size L with 6 
attractive interactions. If we call Xa (for 1 < a < n) the positions of the 
particles, the Hamiltonian of this system is 
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where 7 is the strength of the attractive (7 > 0) interactions. The main goal of 
the present work is to define and to calculate the ground state energy Ef^in, L, 7) 
of (|l|) when n is not an integer (especially when n is small). 

This system of particles in one dimension with 5 interactions has a long his- 
tory in the theory of exactly soluble models |l|, ||, |[ ^, ||. It was first introduced 
to describe a Bose gas by Lieb and Liniger who calculated by Bethe ansatz the 
ground state energy and the excitations for repulsive interactions (that is for 
negative 7) in the thermodynamic limit (n and L go to infinity keeping n/L 
constant) y. 

The problem arose also in the theory of disordered systems: the calculation of 
the free energy of a directed polymer in a random medium in 1 + 1 dimensions by 
the replica method reduces]^, 0, |[ ||, |l^, [ill to finding the ground state energy 
of (|^): if Z{x, t) is the partition function of a directed polymer joining the points 
(0, 0) and (x, t) on a cylinder with periodic boundary conditions {x + L = x) 
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where the random medium is characterised by a Gaussian white noise rjijj^t) 

{v{y,t)viy',t'))^jS{y~y')Sit^t'), (3) 
then the integer moments of Z{x,t) are given for large t byfl^ 



lim — In 

i^oo t 



{Z^\x,t)) 



{z{x,t)y 



-£'o(n,£,7), 



(4) 



where ( ) denotes the average over the random medium and Eo(n, L,"f) is the 
ground state energy of (|^). The knowledge of Eo{n, L, 7) determines for large t 
the whole distribution of In Z(x, t). For example, the variance of In Z{x, t) is 
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Of course, to obtain this variance or other characteristics of the distribution 
of In Z{x,t), one should be able to define and to calculate the ground state 
energy of (Q) not only for integer n, but for any value o/n|l^, |l^ |l^, ^ [l7| . 
Moreover, because of (^), the interactions in (Q) must be attractive (7 > 0); So 
in contrast to the Bose gas initially studied the interactions are attractive 
and the interesting limit is no longer the thermodynamic limit ri — > cx3 but 
rather the limit n — > 0. 

For integer n and L — 00, the n particles form a bound state at energy M, p[ 



Eo{n,oo,'y) = -7' 
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Using this formula for non- integer n helped to understand several properties 
of the distribution of In Z{t) when L is infinite. There are however a number 
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of difficulties with for non-integer n, in particular a problem of convexity: 
c?^ ln(Z")/dn^ should be positive for all n, and so (|[ ||) cannot be vahd at least 
for negative n. We believe that these difficulties are due to the exchange of limit 
t oo and L oo and this is why we try in the present work to determine 
EQ{n, L,j) for finite L. 

The paper is organised as follows: in section 0, we recall the Bethe ansatz 
equations which give the ground state energy of MT) for an (integer) number n 
of particles and we write the integral equation ( |l3[ ) which is a way of solving 
the coupled non-linear equations of the Bethe ansatz. The main advantage of 
this integral equation is that both the strength 7 of the interactions and the 
number n of particles appear as continuous parameters. In section js] we solve 
(^3|) perturbatively in c (where c — 7^/2) for arbitrary n. We notice that the 
coefficients in the small c expansion of Egin, L, 7) are all polynomials in n, thus 
allowing to define the expansion even for non-integer n. In section ^ we show 
how to generate a small n expansion of the solution of (|l3|). We give explicit 
expressions up to order of EQ{n, L,j) and we notice that the coefficients of 
the small n expansion have in general a zero radius of convergence in c. 



2 The Bethe ansatz equations 

The Bethe ansatz consists in looking in the region < a;i < • • • < a;„ < i 
for a ground state wave function ^'(0:1, . . . ,Xn) of the form 

*(xi, . . . , Xn) =Y.^P e*('i^^(i'+-+«"^^(")), (7) 
p 

where the sum in (^ runs over all the permutations _P of {1, . . . , n}. The value of 
^ in other regions can be deduced from (|^) by symmetries. One can show||, |4[ || 
that (|^) is an eigenstate of (|l|) at energy 

i?(n,L,7) = -^ ^ ql, (8) 

1< Q < n 



if the {qa} are solutions of the n coupled equations 

- 9/3 + c 

where 



n ^£i^i£±£, (9) 



(10) 



(A derivation of (p|) can be found in Note that ikj and c i n ||3| become 
here respectively j^qj and —7; the c in |^ and our c defined in (|10[) are thus 
different.) Moreover, for 7 7^ 0, all the qa are distinct. 

There are a priori many solutions of (^. We look for the ground state, 
that is the solution {qa} for which (0) is minimal. When c = 0, the problem 
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reduces to n non-interacting particles {qa} = {0} (we have periodic boundary 
conditions). Because the ground state solution is not degenerate, the solution 
{qa} of (|^) must have the symmetry {qa} ~ {—qa}: depend continuously on c, 
and vanish as c — s- 0. 

Let us introduce the following function of {qa}- 



B(u) = -e 
n 



,<3=(«-l) 



(11) 



where the parameters p{qa) are defined by 

qa 



PiQa) 
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IP 



Qa - qp 
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The function B(u) is a rather complicated (but easier to manipulate than the qa) 
symmetric function of the ground state solution {qa} of (|^). As shown in the 
appendix, it satisfies the integral equation 



/•U 

B{1 + u)- B{1 -u) = nc / e-5("'-™)B(l 
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and the following two conditions 



5(1) = 1, 
B{u) = B{-u). 



v)B(l +u- v)dv, (13) 



(14) 
(15) 



Moreover, the energy (g) can be deduced from of B{u) by 

^2 



Eo{n,L,j) 
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The derivation of (|l3|, |l^, |T5|, |T6|) is given in the appendix. How these 
relations lead to small c or small n expansions is explained in sections ^ and ^. 



3 Expansion in powers of c 

One could try to solve the equation s (p| ) perturbatively in c but the approach 
turns out to be quickly complicated ]19| . Instead we are going to see that the 
integral equation (^|) is very convenient to obtain EQ{n, L,j) for small c. 

It is known[^ (and easy to check from (J^) that the qa scale like -/c for 
small c. Therefore each coefficient Bi{u) of the small c expansion of B{u) 
defined by ( pi] ) is a polynomial in u. 

B{u) = Bo{u) + cBi{u) + c^B2{u) + . . . (17) 

Moreover, conditions ( p^ ) and ( p^ ) impose that all the Biiu) are even, that 
Bq[1) = 1 and Bi[l) = for any i > 1. 
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At zero-th order in c, we find, using (p^): 

Boil + u) - Bo{l - u) ^ 0. (18) 

Thus Bo(u) and Bo{l+u) are both even functions of u. As Bq(u) is a polynomial 
and -Bo(l) = I7 the only solution is 

Bo{u) = 1. (19) 

We put this back into ([l3| ) and get at first order in c 

Bi{l + u) - Bi{\ - u) ^ nu. (20) 

Again, using the fact that Bi{u) is an even polynomial such that -Bi(l) = 0, 
the only possible solution is: 

= ^(w^-l). (21) 

It is easy to see from ([l^ ) that at any order in c, we have to solve 

Bi{l + u) — Bi{l — u) = "some polynomial odd in u" , (22) 

and that there is a unique even polynomial solution satisfying Bi{l) = 0. One 
can generate as many Bi{u) as needed to obtain B{u) up to any desired order 
in c. 

B{u) = l + ^-^ - + — '- + 0{d'). (23) 

Relation (^6|) then gives the energy. Up to the fourth order in c, we find 

^ , N / d c (? n n f m? n \ , \ 

-Eoin, L, ,) . -nin - D 3 + + J^-^^ + (j^ - j + - j ■ 

(24) 

4 Solution for small n 

It is clear from section ^ that if we stop the small c expansion of B{u) at a 
given order, B{u) and Eo{n, L,^) are polynomials in n. This allows to define 
the small c expansion of B{u) or of i?o("-, L, 7) for an arbitrary value of n. 

Moreover, we can write a small n expansion of B{u) by collecting all the 
terms proportional to n'' in the small c expansion of B{u) and calling the series 
obtained bk{u). Then, 

B{u) = l + nbi{u) + n^b2{u) + ... (25) 
Conditions ( p^ , p^ ) impose that bkiu) — bk{—u) and 6^(1) = for all A; > 1. 
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We are now going to describe a procedure which leads to a recursion on the 
bk{u) and allows to write them not only as power series in c but as explicit 
functions of c and u. If we insert ( ]25| ) into (13), we get at first order in n 

6i(1 + m)-6i(1-u) = c / e-it'''-™)^?;. (26) 



It can be checked that a solution of (g6|) compatible with the conditions bi (u) = 
bi{—u) and fei(l) = is 

^ /•+°° cosh ^ -cosh ^ ,2 ^ ^ 

Jo sinh ^ 

There are however other solutions to the difference equation (^): one could 
add to (27) an arbitrary function F(u, c) even and periodic in u of period 2 and 
vanishing at it = 1. If we require that each term in the small c expansion of 
bi{u) is polynomial in u (as justified in section we see that all the terms of 
the small c expansion of F{u, c) must be identically zero. For example F{u^ c) = 
exp(— c~^/'')(cos(7rw) + 1) is an acceptable function. This already shows that 
bi{u) given by ( p7| ) has indeed for small c expansion the series obtained by 
collecting all the terms proportional to n in the small c expansion of section ^. 

If the solution bi{u) (27) had a non-zero radius of convergence in c, it would 
be natural to choose the only bi{u) which is analytic in c at c = by taking 
F{u, c) — 0. UnfortunatcW, this is not the case: by making the change of vari- 
able A^ = /i, expression (|27| ) appears as the Borel sum of a divergent series [po| . 

We found no conclusive reasons why (|2^ ) is the solution of ( p6| ) we should 
select. However, one can notice that, when n is an integer, B{u) is analytic in u 
and goes to zero when u ±ioo (see (|lT|)). Here, the bi{u) given by (^) grows 
like In |u| when u — > ±ioo. Adding a non-zero periodic F{u, c) would either lead 
to an exponential growth in the imaginary direction or introduce singularities 
in the complex u plane. So (|27| ) is the solution of (p6|), analytic in the whole u 
plane, which has the slowest growth in the imaginary direction. 

The same difficulty of selecting the right solution appears at every order in 
the expansion in powers of n. We are now going to explain the procedure we 
have used to select one solution. If we insert ( p5| ) into (|l3|), we have to solve at 
any order k in the small n expansion 

bk{l + u)-bk{l-u) = Mu), (28) 

where (j}k{u) is a function odd in u which can be calculated from the previous 
orders 

4>u{u)=cy^ e'^'^''"-"-^b^{l-v)bk-^-^{l + ^L-v)dv. (29) 

(For consistency, we use 6o(w) — 1.) It can be checked that a solution to ( p8| ) is 

/•+°° cosh cosh ^ , ^ 

hk{u)^ / d\ 2 ^afc(A), (30) 

JQ sinh 
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where afe(A) is given by 



du sm — (pk 



lU 



(31) 



Indeed, the verification is a simple matter of algebra. (We convinced ourselves 
that &fe(ii) ^ \u\ and 4>k{u) ^ In*^""'^ l^^l/"*^ as m ^ ibioo, and that ak{^) ~ 
\n^-^ |A| for A ^ and afc(A) ~ exp(— A2(/c + 1)/4A;) for A — > oo, so that all the 
integrals in (30, ^ converge.) 

As for &i(u), one could add to the &fe(u) given by (|3^, ^l|) an arbitrary func- 
tion Fk{u,c), even and periodic in u of period 2 to obtain the general solution 
of (|2^). However, to be consistent with the small c expansion of section ||, the 
small c expansion of Fk{u,c) should be identically zero. Moreover if we want 
the solution of (|8|) to be analytic in the whole u plane and not to grow too fast 
when u ±ioo, we must take Fk{u, c) — 0. 

Expression for bi{u) is in fact a particular case of the procedure ( P0| , 
^l|) ; when applied to (p6|), it gives indeed ai(A) = ■y/cexp(— A2/2). 

At second order in the small n expansion, we find for A > 



a2(A) 



A 



^2 2 cosh ^ 







iiie"^ - 2 
tanh^ ' A"'^'' tanh^j 



+ 00 

e~ 



(32) 



The expressions of fo3(u) or 03 (A) would be much longer to write and higher 
orders even more complicated. Recursion (^, |o[ ^) allows nevertheless to 
calculate in principle the whole expansion in powers of n. 

Using relation (16) and the expressions ( p7| ) and ( ^ of h\(u) and 02 (A), we 
find that the energy Eo{n, L, 7) is given up to order n'^: 



-Eo{n,L,j) = n - + 



c 

12 



r.3/2 



X3 ^2 
dA -yri 

tanh ^ 



(33) 



+°°A2e- 

aA 

tanh 



2 



^ _i£2cosh4^-2 
d^e 



tanh 



+00 



tanh 



As explained in the introduction, this small n expansion of Eo{n, L, 7) gives the 
cumulants of the free energy in the directed polymer problem. Of course, if we 
expand (p3) in powers of c, we recover (kJ 



5 Conclusion 

In this work, we have developed a method allowing to calculate perturbatively 
the ground state energy of (|l]) for a non-integer number n of particles. We first 
generated for integer n a perturbation series in powers of the interaction c. Each 
term of this series is polynomial in n, allowing to define a small c expansion of 
the energy for non-integer n. This series, at least for small n, has in general 
a zero radius of convergence, in contrast to integer n for which the radius of 
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convergence of the perturbation theory is non-zero pl||. (For n = 2, the closest 
singularities of Eo{2, L, 7) in the c plane lie at c ~ 3.30 ± i 4.12.) 

We believe that the fact that each term in the perturbation theory is poly- 
nomial in n is generic and would be true for an arbitrary pair interaction and in 
any dimension. As the link to directed polymers is valid in any dimension, 
it would be useful and interesting to try to recover our results by doing a direct 
perturbation theory of the Hamiltonian instead of our Bethe ansatz approach 
(which is limited to 1 -I- 1 dimensions and to a (5 potential) in order to see whether 
the calculations could be extended to higher dimensions. 

Our calculation of the ground state energy for non-integer n is based on the 
integral equation ( |l3|) and the conditions (^, |l^). When we tried to solve the 
problem for small n, at each order we had to select a particular solution of a 
difference equation. We did not find a conclusive reason to justify the solution 
we selected, apart from some analyticity properties and growth criterion in the 
complex plane of the variable u. It would certainly be interesting to justify our 
choice (^3|) by calculating the second and the third cumulants of In Z directly 
(and not only perturbatively to all orders in c). 

In our small n expansion of section ^, the terms become quickly very compli- 
cated. There is however a regime, which corresponds to the large c limit of re- 
cursion (H, |0|, 15) where one can handle all orders in the small n expansion pof. 
This allows one to calculate the whole distribution of ln(2'(x, t))/t when t is very 
large and {l/t)ln{Z{x,t)/ {Z{x,t))) of order 1/L. One can then recover [p9[ the 
same large deviation function as found for the asymmetric exclusion process p2|, 
p3| , p4| , ^ , as expected since the directed polymer problem in 1 -f 1 dimen- 
sions and the asymmetric exclusion process are both representatives of the KPZ 
equation p?!, in dimension 1. This strengthens the conjecture that the 

solutions to the difference equations we selected in section ^ give indeed the 
right non-integer moments of the partition function. 

From the point of view of the theory of disordered systems, our results give 
one of the very few examples for which the distribution of Z can be calculated 
exactly. In particular they could provide a good test of the replica approach 
and of other variational methods [Q, ||, |ll|. 

A simple and interesting phenomenon visible in the present work (which is 
generic of all kinds of disordered systems with Gaussian disorder) is that the 
weak disorder expansion (here small c expansion) of non-integer moments of the 
partition function has a zero radius of convergence whereas integer moments 
have a non-zero radius of convergence. This is already visible in the trivial 
example of a single Ising spin a = ±1 in a random Gaussian field h; the partition 
function at temperature T is Z = 2cosh(ft,/T), and it is easy to check that all 
non-integer moments of the partition function have a zero radius of convergence 
in 1/T. 

Acknowledgements We thank Frangois David, Michel Gaudin, Vincent Pasquier, 
Leonid Pastur, Herbert Spohn and Andre Voros for useful discussions. 
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A Derivation of ([I^, [ig, 

Let us first establish some useful properties of the numbers p{qa ) defined by (|l^. 
If the Qa are the n roots of the polynomial P{X) defined as 



P{X) = l[{X-qo,), 



it is easy to see that the p{qa) defined in (|12|) satisfy 
P{X + c) _ ^ p{q^) 

P{X) ^'^X-q^- 



(A.l) 



(A.2) 



(The two sides have the same poles with t he s ame residues and coincide at 
X oo.) Expanding the right hand side of (A.2) for large AT, we get 



P{X 



P{X) 



P(ga) 

X 



X A2 



1 



A4 



(A.3) 



On the other hand, using A.l) and the symmetry {qa} ~ {^qa} we have 

P{X) = A" + ^Eo{n, L, 7)A"-2 ^ o(xn-4)^ (^.4) 

so that 

PiX + c) 



P{X) 
Comparing 



A A2 A3 



1 

A4 



) and ( [A. 51) , we get the relations 

^p{qa) = 

go 

^qapiqa) = , 

^ 2 , ^ 2 /^"-A Eo{n,L,j)L^ 



<3a 



(A.5) 

(A.6) 
(A.7) 

(A.8) 



Moreover, by letting A = — c in 



we get for any root of P{X) 



qa ~ qp + c 



P{(la] 



9/3 



(A.9) 



qa qc. 

Lastly using the symmetry {qa} — {^qa} and the definition (p^), the Bethe 
ansatz equations (^) reduce to 

e'°p(-ga) - e-«°pfe) = 0. (A.IO) 

From the definition (llh of B{u) and the properties ( A.6- |A.l"o| ), it is straight- 
forward to establish (p^ '|16|): the integral equation (|l3) is a direct consequence 
of (|ll|) and ( [All ). Properties (|l|, |^) follow fronT pl, [A^ ) and (|ll|, |a1o| ) 
respectively. Lastly (|l^) is a consequence of A.6 -A.J 
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